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The temporal evaluation of the 27-day variation of the three dimensional (3D) galactic cosmic ray 
(GCR) anisotropy has been studied for 1965-2014. 3D anisotropy vector was obtained based on 
the neutron monitors and Nagoya muon telescopes data. We analyze the 27-day variation of the 
(1) two dimensional (2D) GCR anisotropy in the ecliptic plane, and (2) north-south anisotropy 
normal to the ecliptic plane. Studying the time line of the 27-day variation of the 2D GCR 
anisotropy, we confirm that the average amplitude in the minimum epoch of solar activity is po¬ 
larity dependent, as it is expected from the drift theory. The amplitude in the negative polarity 
epochs is less as we had shown before. The feeble 11-year variation connected with solar cycle 
and strong 22-year pattern connected with solar magnetic cycle is visible in the 27-day variation 
of the 2D anisotropy for 1965-2014. We show that the 27-day variation of the GG index (being 
a measure of the north-south asymmetry) varies in accordance to solar cycle with a period of 
11-years, being in good agreement with the 27-day variation of the At component of the GCR 
anisotropy calculated by the IZMIRAN group. Detailed analysis are presented for the solar mini¬ 
mum 2007-2008 of solar cycle no.23 and solar maximum 2013-2015 of solar cycle no 24. In the 
solar cycle no. 24 GG index, calculated by Nagoya telescopes data, is highly anticorrelated with 
By component of the interplanetary magnetic field (IMF) and shows a clear recurrent changes 
related to the Sun’s rotation. 
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1. Introduction 

An average solar diurnal anisotropy of galactic cosmic rays (GCRs) is formed due to modula¬ 
tion of GCR particles by the solar wind. It has been explained based on the diffusion-convection 
theory of GCR propagation in the heliosphere [P-[^, as a consequence of the equilibrium estab¬ 
lished between the radial convection of GCR particles by solar wind and the inward diffusion of 
GCR particles (owing to an existence of radial gradient) along the interplanetary magnetic held 
(IMF). Despite, three-dimensional (3D) anisotropy vector determines the distribution of the stream 
of cosmic rays in the 3D heliosphere, nevertheless, the complete nature of the 3D GCR anisotropy 
from top to bottom is not understood yet. Scientihc staff of IZMIRAN’s cosmic ray laboratory 
(http://hehos.izmiran.troitsk.ru/cosray/main.htm), have calculated the components Ar, Af, and At 
of the 3D anisotropy by Global Spectrographic Method (GSM) 0, [||] based on hourly data from 
all operating neutron monitors (NMs). Unfortunately, the derivation of the At component is pos¬ 
sible with an accuracy up to constant, and so, a value of the north-south anisotropy At obtained 
by GSM method is not accurate [http : / /cr20.izmiran.rssi.ru/AnisotropyCR/index.php\. That is, 
sorrowfully, a deficiency of the GSM related to the nature of NMs data. In general determination 
of the two dimensional (2D) anisotropy in the ecliptic plane is feasible based on establishing the 
radial Ar and tangential Af components by the harmonic analysis method for an individual detector 
(e.g., NM or Muon Telescope (MT)). Alania et ah, show that GSM gives the same results as the 
individual NM with cut-off rigidity < 5 GV. 

Swinson [Q] hrst suggested that the north-south anisotropy could be related to the cosmic 
ray flow caused by a positive heliocentric radial density gradient Gr of cosmic rays and the By 
component of the interplanetary magnetic field (IMF) B, as expressed by the vector product. By x 
Gr. Since the average By and Gr he in the ecliptic plane, the direction of the vector product By x Gr 
(or the north-south anisotropy) is expected to be perpendicular to the ecliptic plane and to reverse 
direction with changes of the magnetic field direction. It is directed upward when By is positive 
(IMF away from the Sun-A > 0) and vice versa when By is negative (IMF toward the Sun-A < 0). 
Papers [^- [^ studied the cosmic ray anisotropy vector in three dimensions using data from polar 
located NMs. 

An alternative method to study north-south anisotropy is proposed by Mori and Nagashima 
I®. They have introduced an index GG calculated from Nagoya MT data as: GG = \ [{A9N — 
495) + (49A — 49£')]. The terms represent the counting rate differences for telescopes pointing in 
North, South, and East directions at 49° zenith angle. The GG index is free of noise in isotropic 
intensity caused by Forbush decreases, periodic variations, atmospheric temperature effects, and 
geomagnetic cut-offs. GG index mainly reflects change of difference between the intensity from 
north polar direction and that from north and parallel to the equatorial plane directions. Counting 
rate differences of the north-south and north-equatorial could not be precisely contained the same 
type of information. In spite, a GG is accepted by cosmic ray community as a good alternative 
index to study the north-south asymmetry of the cosmic ray flux e.g., dnii- 

The large amount of papers dedicated to the recurrent variations of cosmic rays related to the 
period of the solar rotation, further in this paper called 27-day variation, are devoted to the 27-day 
variation of the GCR intensity. The 27-day variation of the GCR anisotropy in general has been 
studied less intensively up to present. It partially is connected with the small amplitudes of the GCR 
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Figure 1: Time line of the amplitudes of the 27-day variation of the 2D GCR anisotropy (A27A), smoothed 
over 13 Sun’s rotations, for Climax, DeepRiver, Kiel, Moscow and Oulu individual NMs and for all com¬ 
bining NMs by GSM method for 1965-2014. 


anisotropy (< 0.3%, measured by NMs) and with a large dispersion comparable with the accuracy 
of hourly data of NMs. The existence of the 27-day recurrence of the GCR anisotropy was first 
reported by Yoshida and Kondo [1^] and later confirmed by ofhers. The 27-day variation of fhe 
norfh-soufh anisofropy was sfudied in fhe series of papers by Swinson and coaufhors [13|-[p^]. 
They showed fhaf 27-day variafion of fhe norfh-soufh anisofropy is correlafed wifh solar acfivify 
and fhis correlafion is nof clearly dependenf upon solar magnetic polarity. Recenfly, Yeeram el 
al., [15] sfudied recurrenl Irains of enhanced GCR anisofropy under influence of corofafing solar 
wind sfrucfures near fhe recenf solar minimum. They have nofed fhaf during fhe recurrenl lime 
(~ 27 days) fhe GCR anisofropy was often enhanced or suppressed in accord with the sign of the 
interplanetary magnetic field B, possibly demonslraling a conlribulion from a mechanism involving 
a soufhward gradienf Gt in fhe GCR densily, and an anisofropy caused by driff, BxGt. However, 
fo be sure in fhis suggestion Ihere is necessary fo carry oul an additional quanfitafively analyze 
involving dala of NMs wifh low cul-off rigidilies (< 5 GV). 

Our aim in fhis paper is fo sludy long-lerm changes of fhe 27-day variation of fhe GCR 
anisofropy for 1965-2014. Defailed analysis are presented for the solar minimum 2007-2008 of 
solar cycle no.23 and for the solar maximum 2013-2015 of solar cycle no. 24. 


2. Long term changes of the 27-day variation of the GCR anisotropy 

Generally the vector of the GCR solar anisotropy A, describing the stream of cosmic rays in the 
heliosphere, is three dimensional (3D). In this part of paper we investigate the two dimensional (2D) 
case of the GCR anisotropy vector being the projection of the vector A on the ecliptic plane (Ar, 
A/). Using harmonic analysis method we calculated the amplitudes of the 2D 27-day variations 
of the GCR anisotropy {AllA) for Climax, DeepRiver, Kiel, Moscow and Oulu NMs for each 
Sun’s rotation period during 1965-2014. Amplitudes A27A of the 2D GCR anisotropy for individual 
NMs are compared with the 27-day variation of the GCR anisotropy vector (for 10 GV particles) 
calculated using Global Spectropgraphic Method (GSM) by the cosmic ray laboratory of IZMIRAN 
combining all NMs (http://helios.izmiran.troitsk.ru/cosray/main.htm). Fig. |I| presents for the first 
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Figure 2: Daily GG index. By component and the sector structure of the IMF in 2008. 



Figure 3: Temporal changes of the smoothed over 13 Sun’s rotations amplitudes of the 27-day variation of 
the GG index {A21GG, solid line) and At component of the 3D anisotropy (A21At, dashed line). A21GG is 
presented for 1970-2014, A27Af - for 1970-2006. 


time the temporal changes of the 27-day variation of the 2D GCR anisotropy A27A smoothed 
over 13 Sun’s rotations for 1965-2014 for Climax, DeepRiver, Kiel, Moscow and Oulu individual 
NMs and for all combining NMs by GSM. We confirm our previous finding [^] that near the solar 
minimum the amplitudes of the 27-day variations of the 2D GCR anisotropy determined using 
GSM basically do not differ from the amplitudes found by harmonic analysis for an individual NM 
with cut-off rigidity < 5 GV. Figure |I] shows that this feature is kept for the whole solar cycle and 
solar magnetic cycle. The time lines of the 27-day variation of the 2D GCR anisotropy confirm fhaf 
fhe average amplifude in fhe minimum epoch of solar acfivify is polarity dependenf, as if is expecfed 
from fhe driff fheory. The amplifude in fhe negative polarify epochs is less as we had shown before 
1^]. The feeble 11-year vaiiafion connecfed wifh solar cycle and sfrong 22-year pattern connected 
wifh solar magnefic cycle is visible in fhe 27-day variation of fhe 2D anisofropy for 1965-2014. 

The nexf goal is fo sfudy fhe periodic characfer of fhe norfh-soulh componenf of fhe GCR 
anisofropy. If is performed by fwo differenf approaches: (1) using At componenf obfained by 
GSM and (2) using GG index being fhe measure of fhe norfh-soufh anisofropy. Subsequenfly using 
harmonic analysis mefhod we calculated the amplitudes of the 27-day variations of the GG index 
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{AllGG) and Af component {AllAt) being the measure of the north-south anisotropy normal to the 
ecliptic plane. We compare results for the GG index with analysis of the At component obtained 
by GSM. As far the formulation of the north-south anisotropy is based on the drift model, the GG 
index being the measure of the north-south asymmetry is highly anticorrelated with By component 
of the IMF. As an example Fig. ^ presents daily GG index, By component and the sector structure 
of the IMF in 2008. Figure ^ presents time lines of the amplitudes of the 27-day variations of the 
GG index and the At component of the 3D GCR anisotropy for 1970-2014; AllGG is presented 
for 1970-2014, All At - for 1970-2006 (data of the GG index starts in 1970 and At component 
after 2006 is not available at IZMIRAN website). The general behavior of the 27-day recurrence 
of the GG index and normal component At of the 3D anisotropy is similar. We show that the 27- 
day variation of the GG index varies in accordance to solar cycle with a period of 11-years, being 
in good agreement with the 27-day variation of the At component of the GCR anisotropy. Fig. 
^ shows that there is a clear positive correlation between the changes of the AllGG and All At, 
indicating with a high probability about the existence of theirs common source. Thus, data from 
cosmic ray detectors at low (NM) and high (MT) rigidities confirm no significant dependence of 
the 27-day variation of the north-south anisotropy upon solar magnetic polarity. 

3. 27-day variations of the GCR intensity and 2D anisotropy, GG index, solar wind 
velocity and IMF components in the Solar Cycle no. 24 

The current Solar Cycle No. 24 (SC 24) is the least active in a Space Age and even in the 
perspective of a century. The period, from 2007 up to now, because of its very peculiar nature, is 
broadly investigated in many aspects. One of them is the long-lasting very deep solar minimum 
2007-2009 with well established recurrent structures related to the Sun’s rotation period (~27 
days). Recently, the peak of the SC 24 has arrived and the second-half of the solar cycle is officially 
sfarfed. Also in fhe maximum of solar acfivify in 2014-2015, fhe 27-day recurrenf variations were 
clearly visible in cosmic ray infensifies, GG-index, solar wind velocity and heliospheric magnefic 
field. So fhe SC 24 offers fhe possibility fo sfudy 27-days recurrenf variafions over a wide range of 
condifions: from low fo high solar activity. 

To reveal fhe quasi-periodic effecfs connecfed wifh fhe rofafion of fhe Sun in fhe SCs 23 and 24, 
we employ fhe wavelef fime-frequency specfrum fechnique developed by Torrence and Compo 
In our calculation we used fhe Morlef wavelef mofher function. We performed wavelef analysis for 
fhe daily solar wind velocify V, By componenf of fhe IMF, GCR infensify for Oulu NM, GG index, 
Ar and Af componenfs of fhe 2D GCR anisofropy during 2007-2008 solar minimum and 2013- 
2015 solar maximum, considering as a sampling lime interval equaling fwo years. In our case a 
lime inlerval of fwo years gives a good enough slalislics as far our aim is fo reveal recurrences < 
27-day period, which is ~ 4% from whole sampling period of 26 solar rolalions. 

Resulfs of calculation using wavelef analysis melhod for V, By and GCR intensity are pre- 
senled in Fig. ^ (solar minimum) and Fig. ^ (solar maximum). Resulfs of analysis for Ar and 
Af componenfs of fhe 2D GCR anisofropy for Oulu NM and GG index are presented in Fig. || 
(solar minimum) and Fig. ^ (solar maximum). Figures ^ and ^ panel (c) presenl very clear 
quasi-periodic changes in GG index relaled fo fhe Sun’s relation (~ 27 days) for almost whole 
analyzed time intervals. Similar quasi periodic character is clearly visible in V and By component 
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Figure 4: Wavelet analysis of the daily solar wind velocity V (a), By component of the IMF (b), GCR 
intensity for Oulu NM (c) for 2007-2008. 



Figure 5: Wavelet analysis of the Ar (a) and Af (b) components of the 2D GCR anisotropy and daily GG 
index (c) for 2007-2008. Dashed line designates the period of 27 days. 


of the IMF. 27-day variation of GCR intensity is also well established at the end of 2007 and in 
2008 during solar minimum (Fig. ^ and at the end of 2014 and in the beginning of 2015 during 
solar maximum (Fig. |^. Although recurrent variations connected with corotating structures (~ 27 
days) are clearly established in almost all solar wind parameters, 2D GCR anisotropy shows a weak 
27-day variation, only in some periods. This is connected with large dispersion of daily Ar and Af 
components of the 2D GCR anisotropy. 

4. Conclusions 

1. The time line of the 27-day variation of the 2D GCR anisotropy confirm that the average 
amplitude in the minimum epoch of solar activity is polarity dependent, as it is expected 
from the drift theory. The amplitude in the negative polarity epochs is less as we had shown 
before. The feeble 11-year variation connected with solar cycle and strong 22-year pattern 
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Figure 6: Wavelet analysis of the daily solar wind velocity V (a), By component of the IMF (b), GCR 
intensity for Oulu NM (c) for 2013-2015. 



Time (days) 



Figure 7: Wavelet analysis of the Ar (a) and Af (b) components of the 2D GCR anisotropy and daily GG 
index (c) for 2013-2015. Dashed line designates the period of 27 days. 


connected with solar magnetic cycle is visible in the 27-day variation of the 2D anisotropy 
for 1965-2014. 

2. We show that the 27-day variation of the GG index varies in accordance to solar cycle with 
a period of 11-years and is almost independent of solar magnetic polarity, being in good 
correlation with the 27-day variation of the At component of the GCR anisotropy calculated 
using the GSM method. 

3. We study recurrent variations in 3D GCR anisotropy connected with corotating structures 
observed in the heliosphere in the Solar Cycles no. 23 and 24. 2D GCR anisotropy generally 
does not show evident 27-day variation, but in some periods weak recurrent character is 
visible. This is connected with large dispersion of daily Ar and Af components of the 2D 
GCR anisotropy. Using wavelet time-frequency method we reveal clear 27-day waves in 
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the GG index, being the measure of the north-south anisotropy, for almost whole analyzed 
period. 
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